Mon. Not. R. Astron. Soc. 000,[T]-xxx (2010) Printed 20 January 2013 (MN BTfiX style file v2.2) 



The Structure of Molecular Clouds: III - A link between cloud 
structure and star formation mode 

Jonathan Rowles x * Dirk Froebrich : f 

1 Centre for Astrophysics & Planetary Science, The University of Kent, Canterbury, Kent CT2 7NH, U.K. 



Accepted Received ; in original form 



ABSTRACT 

We analyse extinction maps of nearby Giant Molecular Clouds to forge a link between 
driving processes of turbulence and modes of star formation. Our investigation focuses on 
cloud structure in the column density range above the self shielding threshold of 1 mag Ay and 
below the star formation threshold - the regime in which turbulence is expected to dominate. 

We identify clouds with shallow mass distributions as cluster forming. Clouds that form 
stars in a less clustered or isolated mode show a steeper mass distribution. Structure functions 
prove inadequate to distinguish between clouds of different star formation mode. They may, 
however, suggest that the turbulence in the average cloud is governed by solenoidal forcing. 
The same is found using the A-variance analysis which also indicates that clouds with a 
clustered mode of star formation show an enhanced component of compressive driving in the 
turbulent field. Thus, while star formation occurs in each cloud, independent of the turbulent 
driving mechanism, compressive forcing appears to be associated with the formation of stellar 
clusters. 
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1 INTRODUCTION 

Stars form due to the local collapse of material in molecular clouds. 
The conditions prior to the collapse result from a complex inter- 
play between self-gravity, turbulence, magnetic fields and thermo- 
dynamics. Understanding the effects of each of these influences 
leads to a better knowledge of how stars form. In this paper we 
are interested in the role played by interstellar turbulence and it's 
effect on the structure of Giant Molecular Clouds (GMCs). 

The subject of astrophysical turbulence is complex and not 
fully understood, though there are many useful reviews on the topic 
(e.g. Elmegreen & Scalo ( 12004b . Scalo & Elmegreen (20041, Bran- 
denburg & Nordlund (2009 ))- Here we are interested in comparing 
the structure of nearby GMCs to investigate whether it has an in- 
fluence on the observed star formation properties. Ultimately this 
might also be used to determine the nature of the turbulent field, i.e. 
driven by a compressive forcing or solenoidal driving (e.g. Feder- 
rath et al. dMot ). 

We will probe the column density structure in nearby GMCs 
by means of extinction maps derived from near-infrared obser- 
vations. This is the least biased way to estimate column density 
(Goodman et al. ([2009)). In the first paper in this series (Rowles 
& Froebrich (2009), hereafter Paper I), we presented new all-sky 
extinction maps derived using data from the 2 Micron All-Sky Sur- 
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vey (2MASS, Skrutskie et al. J2006b ). We used the median near- 
infrared colour excess technique (NICE) to calculate the extinction 
(see Lada et al. ( 11994b ). The nearest 25, 49 and 100 stars to the 
centre of each pixel were used, hence the noise can be considered 
constant throughout the map. These extinction maps are therefore 
referred to as con-noise maps. 

In our second paper (Froebrich & Rowles (2010), hereafter 
Paper II), we analysed the column density and mass distributions 
of a selection of 16 nearby GMCs. To facilitate this we determined 
new extinction maps using only the stars within each pixel (i.e. no 
oversampling), which therefore have a constant spatial resolution 
and are referred to as con-res maps. As a result of the analysis we 
found a universal star formation threshold of about 6.0 ± 1.5 mag 
Ay. This threshold separates two different regions in the clouds. 
Below the threshold, at low column densities, turbulence dominates 
the structure, while at higher column densities gravity is the dom- 
inant force. The low Ay part of the clouds could be fitted by a 
log-normal distribution. There were significant differences in the 
slopes of the column density and mass distributions when consid- 
ering only the low Ay regions. This shows that the properties of 
the turbulence differ depending on the environment of the cloud. 
Regarding the high Ay regions, we found no such differences, im- 
plying that gravity solely dominates these parts. 

Using our extinction maps we can derive column density 
structure functions similar to velocity structure functions (e.g. 
Padoan et al. ( 2002 )) for each molecular cloud. This will allow us to 
perform a comparison with models of interstellar turbulence. Pre- 
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dictions of structure functions resulting from a turbulent medium 
have been presented e.g. by Schmidt et al. ( 12008b . Kolmogorov 
dT94lT >. She & Leveque (1953} and Boldyrev J2002t (hereafter S08, 
K41, SL94 and B02, respectively). 

In this paper we test for correlations between structure func- 
tion parameters and the properties of the clouds. We also examine 
the cloud structures using the A-variance technique (see Stutzki et 
al. ([19981 and Ossenkopf et al. (2008a I) and in particular the mass 
spectral index scaling coefficient. In Sec.|2]we describe the methods 
used for our analysis. In Sec.[3]we give the results for the clouds se- 
lected. We discuss these results and give conclusions in Secs.|4]and 
[5] respectively. 



2 METHOD 

2.1 Structure functions 

The general definition of the structure function is given in Eq.[T] 
(e.g. Lombardi et al. 1 2008 1; Padoan et al. {2002 1 and Padoan et al. 
(2003 )). Here the equation is expressed in terms of our observable 
- the column density or optical extinction Av- 

S P (Ar) = (\A V (/) - A V (r' - Ar) | p ) (1) 

Ar is the distance between points, r' represents a position in 
the map, Av is the optical extinction at r' (or r' — Ar) and p is 
the order of the structure function. For p = 2 the equation is the 
two-point correlation function of the extinction map. The brackets 
() denote that the average over all pixel positions r' and all possible 
directions for the separation Ar of points is applied. For each order 
p we find the scaling exponent s(p) by fitting a power-law to values 
of Sp(Ar) against Ar. This assumes that the scaling exponents 
are related to the structure functions by Eq.|2] (e.g. Padoan et al. 
( 120031 ): 

S p (Ar) oc Ar s(p) . (2) 

Equation[2] allows us to determine the scaling exponents s(p) 
which are then normalised to the third order s(3), as a universal be- 
haviour should be exhibited at low Reynolds numbers (determined 
by Benzi et al. ( f!993l ). 

For each GMC investigated we used the range of spatial scales 
Ar from 0.1 pc to 1.0 pc to fit the power law exponent, in order 
to be able to compare the results for all clouds. We then follow 
Padoan et al. (2002) and assume the column density scaling ex- 
ponents s(p)/s(3) are equivalent to the velocity scaling exponents 
C(p)/C(3) (Dubrulle (1994)). They can then be expressed by a re- 
lation of the form as shown in Eq.[3] 

Here /3 is the intermittency and A is related to the co- 
dimension C and intermittency by: A = C ■ {1 — ft). The frac- 
tal dimension D of the cloud is related to the co-dimension by 
D — 3 — C. Using this equation, the K41 relation can be expressed 
by setting A = 0. As mentioned previously we use the projected 
column density scaling parameter s instead of the velocity scaling 
parameter £ to keep the same description of the structure function 
as in Padoan et al. (2002), though it is not known a priori whether 
the parameters correspond to each other. However, Padoan et al. 



( 2002 1 found that the scaling exponents derived from integrated in- 
tensity images follow the velocity scaling using the B02 relation. 

We may then derive the parameters A, C, f3 and D for each 
cloud and compare with the velocity field structure functions ob- 
tained from the aforementioned theoretical works by K41, SL94, 
B02 and S08. When fitting our data to Eq.[3] we varied A in the 
range 0.02 to 1.20 with increments of 0.01, and C ranged from 
0.02 to 3.00 with increments of 0.01 to find the best match. Only 
fits with an rms better than 0.1 (in units of s(p)/s(3)) were con- 
sidered a good fit. This is justified given that rms is not an absolute 
measure of the goodness of fit. 

Regarding the previously published models, the simplest to 
consider is that of K41. The energy over a wide range of lengths 
(known as the 'inertial range') in turbulent flows is redistributed 
from larger scales into ever smaller scales until the effects of viscos- 
ity become important. In K41, Kolmogorov considered the struc- 
ture velocity functions S2 and S3, showing that both are power laws 
of the form of Eq.|2] with the exponents being 2/3 and 1, respec- 
tively. Due to the self-similarity of turbulence at different scales, 
this law can be extended to all powers of p such that s(p) = p/3. 

However, experimental measurements for the scaling expo- 
nents show a deviation from the K41 relation for turbulence, for 
higher orders of p - known as intermittency. This is exhibited in the 
non-Gaussian tails when plotting the probability density functions 
(PDFs) of e.g. the column density, and anomalous scaling of higher 
order structure functions (e.g. Anselmet et al. ([T984)). Therefore, 
we need a relation for the scaling that considers the effects of inter- 
mittency. 

SL94 derived a scaling relation that accounts for intermittency. 
Their relation shows good agreement when compared to simula- 
tions of structures that are not influenced by magnetic fields. For 
incompressible turbulence SL94 derived that C = 2 and A = 2/3. 
Therefore, the SL94 model assumes that D = 1, i.e. that the most 
intermittent structures are filaments. 

B02 extended the SL94 relation to model highly supersonic 
turbulence and also take account of magnetic fields. They found 
that C = 1 with A = 2/3 as with the SL94 model. Under the B02 
model it follows that D — 2, i.e. that the structures are sheet-like 
in form. 

In S08, the scaling exponents were found to be well described 
by log-Poisson models. In these models A ~ 1, rather than 2/3. 
The modelling applied two types of forcing to drive the turbulence, 
solenoidal (divergence-free) and compressive (rotation-free). Val- 
ues for the parameters A and C differ depending on the type of 
forcing applied. Perhaps the nature of forcing in real clouds can be 
found by comparing the observational data with the S08 model. 

2.2 The effect of noise on the structure functions 

To ascertain the effect of noise in the data when calculating struc- 
ture functions from an image we undertook a number of tests. In 
particular we desired to know whether (increased) random noise 
leads to a significant and/or systematic deviation of s(p) /s(3) com- 
pared to a test image without noise, or a low noise Ay map. In 
the first instance we made an artificial Av map of 400 x 400 pix- 
els containing a circle of radius 124 pixels, where the pixel values 
in the circle are varied linearly from 0.81 mag (edge) to 10.0 mag 
(centre). This may be considered analogous to a perfectly spher- 
ical cloud with a smooth density variation from the core to the 
edge (admittedly, an unlikely scenario in practice). All other pix- 
els in the map were set to zero. The structure function of this map 
happened to correspond to Kolmogorov type turbulence. We then 
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added 0.28 mag Gaussian noise (corresponding to the 1 a noise in 
our 49*' 1 nearest neighbour map - see Paper I) to this map and fit the 
parameters A, C, f) and D. This process was continued by adding 
further amounts of noise to the basic map with 0.28 mag (1 a) noise. 
The quantities added were a further 0.0625a, 0. 125a, 0. 1875a and 
0.25a, generating five images in total. The results of this exercise 
are discussed in Sec. 13.21 

In addition to investigating the effect of noise on artificial 
clouds, we added random noise to a selection of real clouds. 
The clouds chosen were Chameleon, Circinus, Corona Australis, 
Orion A and Orion B. We used the nearest 49 th nearest neighbour 
Ay maps and added 0.0625a, 0.125a, 0.1875a and 0.25a noise. 
We performed the structure function analysis as with the artificial 
clouds. These results are also discussed in Sec.13.21 



2.3 The A-variance 

The A-variance method for analysing molecular cloud structures 
in astronomical images (i.e. in two dimensions) was introduced by 
Stutzki et al. J1998b and improved by Ossenkopf et al. ( 2008b}. The 
method works by measuring the quantity of structure on a particular 
length scale (e.g. Ar) and filtering the data /(Ar) with a spheri- 
cally symmetric 'down-up-down' type function. The function (de- 
noted Q ; (Ar)), which treats different regimes of Ar separately, is 
given by Eq.[4]below. 

Q(Ar) = 0, (Ar)-Q, (Ar) (4) 

— I > — ^ l, core > — ^ l,ann 

The RHS of Eq.|4]is defined using Eqs.[5]and[6] 



r-\ /a s 4 / A 



(5) 



Q(Ar) 



nl 2 (v 2 - 1) 



exp 



exp 



/ Ar 2 



V(*/2) 2 



The A-variance, o\ ( Ar) is then defined in Eq.[7] 



a A 2 (I) 



/(Ar)*0,(Ar; 



(6) 



(7) 



The filter function (J_) ; (Ar) denotes the Fourier transform of 
the filter function with size / and diameter ratio v. Two particu- 
lar filter types were found to best probe the structure in molecular 
clouds, a 'Mexican hat' filter and a 'French hat' filter (Ossenkopf 
et al. (2008a I). Both filters treat the core and the annulus separately 
(Eq.|4j. The optimal filter function for molecular clouds was found 
to be either a Mexican hat filter with a diameter ratio of 1.5, (or 
a French hat filter with a diameter ratio of 2.3 - Ossenkopf et al. 
(2008a)). We used a Mexican hat filter with a diameter ratio of 1.5 
for our analysis. 

We applied this method to our con-noise maps. To further in- 
vestigate the dimensions of structures in the clouds (i.e. filaments 
etc.) we also applied this method to the corresponding star density 
maps. The structure shown in the star density maps is enhanced 
since the extinction in the Av maps is roughly proportional to the 
log of the star density. Therefore, both the Av and star density 
maps were used for this analysis. The results are shown in Sec. |3.3| 

Using the A-variance technique allows us to also calculate the 
mass spectral index scaling exponent (denoted a) of each cloud. 



This is found by fitting a power-law to a\ against scale. The power 
law is fitted over the range of size of the cloud, i.e. between the 
smallest and largest scales of the cloud. The smallest scale is deter- 
mined by eliminating the pixels affected by oversampling. For the 
largest scale we chose the point where the increase in a\ begins to 
tail off before the peak. These results are also shown in Sec.|3.3| 



2.4 The effect of noise on the A-variance 

As with the structure functions it is also necessary to consider the 
effects of noise on the A-variance results. This was investigated 
in a similar way as for the structure functions (see Sec. |2.2[ (. The 
same set of test images of 400 x 400 pixels containing a circle, 
as well as real clouds were used. The A-variance of these image 
was calculated. Increments of Gaussian noise were added to each 
of the images in exactly the same way as for the structure function 
analysis. The results are presented in Sect.|3.4| 



2.5 Molecular cloud density 

We estimate the average physical density of the cloud material in 
order to relate it to the cloud properties obtained from our analysis 
of the column density and the star formation properties. The aver- 
age density of material can be estimated from the mass of the cloud 
M above a certain Av threshold, the number of pixels N the cloud 
is covering in our Av map with values above this threshold, and 
the distance d to the cloud in the following way: 



p[M Q /pc 3 ] = 



8.78 ■ 10 15 ■ M[M e ] 
Ni ■ (a;["Mpc]) 3 



(8) 



Here the average density is expressed in solar masses per cubic 
parsec, the mass in solar masses, the pixel size x in arcsec and the 
distance d in parsec. The power of 3/2 at the number of pixels of 
the cloud assumes that the two visible dimensions in the plane of 
the sky are a good representation of the dimension along the line 
of sight. Note that if one wants to express the density in units of 
particles per cubic centimeter, then the constant 8.78 • 10 lj needs 
to be multiplied by a factor of 40. 

We determine two average densities for the purpose of this 
paper, i) The average density of the entire cloud. This includes all 
material which is above the self shielding column density threshold 
of 1 mag Av- ii) The average density of the material above the 
cloud's star formation threshold. This corresponds to the density 
of the material which is potentially involved in star formation. The 
results are shown in Sect. 13.51 



3 RESULTS 

3.1 Structure function results 

For the analysis of the structure functions we used the three avail- 
able con-noise maps (see Paper I) for the following reasons: 

i) The noise has a less predictable effect on the results when 
using the con-res maps. This is caused by the fact that each pixel 
has a different number of stars contributing to the Av value, and 
hence variable noise. In some cases only very few stars are used. 
One drawback of using the con-noise maps for this analysis is the 
variable spatial resolution between pixels. This is most notable in 
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Figure 1. A plot showing the variation of s(p) /s(3) against p for the Au- 
riga 1 cloud and the 49 th nearest neighbour map. Our data (+ signs) is com- 
pared against the models of K41 (solid line), SL94 (dotted line) and B02 
(dashed line). The best fit values of A, C, f), D and the rms value are 
shown in the upper left corner. A plot of the difference of the fit from the 
data against p is shown as a red solid line in the lower panel. 



regions of high extinction. However, we only include pixel separa- 
tions above 0. 1 pc (see above), which are generally above the spa- 
tial resolution. Furthermore, we determine the structure functions 
only for the parts of the cloud where turbulence dominates, i.e. 
where Av is greater than 1 mag and less than the cloud's individual 
star formation threshold (generally below 8 mag, see Paper II). 

ii) The obtained values for s(2) and s(3) are very low when 
using the con-res maps. This could be a result of the generally 
higher and, more importantly, variable noise in these maps. With 
the con-noise maps we found much more reasonable values for 
s(2) and s(3). We calculate a mean value for s(2) of 0.8±0.4 and 
s(3) of 1.1±0.6. For individual clouds, using Taurus as an exam- 
ple, we find s(2)=0.98 and s(3)=1.41, compared to s(2)=0.77 and 
s(3)=1.10 (Padoan et al. ( 125531 ). For Perseus we find s(2)=1.24 
and s(3)=1.75, compared to s(2)=0.83 and s(3)=1.18. 

Thus, we determine the structure function of every GMC for 
each of the three available Av maps, utilising only extinction val- 
ues above the self-shielding limit (1 mag) and below the star for- 
mation threshold of the particular cloud. Note that every individual 
extinction value used for these calculations has a signal-to-noise 
ratio of at least three (Paper I). Hence, the obtained structure func- 
tions are highly reliable, as indicated by the small uncertainties for 
s(p)/s(3) (e.g. Fig.[TJ. 

For each of the three structure functions we determine the best 
fitting set of parameters A, C, /3 and D, and a scatter for each 
parameter by considering only the fits with an rms of less than 
0.1 (in units of s(p)/s(3)). The final structure function parameters 
for each cloud are then determined as weighted averages of the 
three obtained values. An example plot showing s(p)/s(3) against 
p for the Auriga 1 cloud together with the best fit is shown in Fig.[T] 
where we used the 49 th nearest neighbour map. Similar plots for 
all clouds are shown in the Appendix |Al| 

The determined parameters of the structure function (A, C, /3 
and D) for all the clouds and their uncertainties are listed in Ta- 
bleQ] The rms of the fit is also shown. The values for A lie in the 
range 0.45 to 1.08 with a mean of 0.85±0. 16. Therefore, even after 
allowing for uncertainties, the values are more or less in the middle 



of those quoted in the literature (A=2/3 in e.g. Padoan et al. (2003 ), 
SL94, B02; A=l in e.g. Schmidt et al. fMM ). 

The co-dimension values C lie between 1.18 and 2.51 with 
a mean of 1.92±0.37. In turn the fractal dimension of our sample 
of clouds has a value of 1.08±0.37, rendering them filament like 
rather than sheet like. 

The values of j3 (the degree of non-intermittency) in our sam- 
ple range from 0.38 to 0.65 with a mean of 0.53±0.08. The litera- 
ture usually quotes a values for f3 of either 2/3 (SL94) or 1/3 (B02). 
Therefore, the value for our GMC sample lies about halfway be- 
tween these values. 

3.2 The effect of noise on the structure functions 

We added noise to a test image and to a selection of the extinction 
maps of real clouds to investigate the effect on the structure func- 
tions and their fitted parameters (A, C, /3, D; see Sec. |2.2} . The 
results of this exercise are shown in Table|2] 

In the top part of the table we list the values of the structure 
function parameters obtained for our test image and the test im- 
age plus noise. As one can see, the artificial cloud corresponds to a 
structure in agreement with Kolmogorov type turbulence (A ~ 0). 
This also implies that the actual values of the other parameters are 
meaningless (see Eq.[3](. Once our typical noise of 0.28 mag A v has 
been added the parameters change significantly. The values for A 
are 0.5 and the structure corresponds to an almost two-dimensional 
object (75=1.9). When the noise is further increased no systematic 
and/or significant changes are observed. This has two implications: 
i) Our typical noise can change the structure function parameters 
significantly compared to a column density map which is free of 
noise, ii) Once the image has some observational noise, a small in- 
crease in the noise will not change the values of the determined 
structure function parameters significantly or systematically. All 
changes are well below the uncertainties of the parameters. This 
implies that it does not matter which of our Av maps we use (and 
indeed justifies that we average the results from all three maps), as 
the variation of a factor of two in the number of stars leads only to 
an increase/decrease of the noise in the maps by 33 % (caused in 
part by the covariance, see Paper I). 

The bottom part of Table|2] illustrates the latter part of the 
above discussion for a real cloud (Corona Australis). With the stan- 
dard noise in the 49 th nearest neighbour map of 0.28 mag Av we 
obtain the listed values for the structure function parameters. When 
increasing the noise by another small amount (at maximum 25 %), 
again we do not find any significant and/or systematic changes in 
the determined structure function parameters. 

3.3 A-variance results 

We applied the A-variance technique (as described in Sec. |2.3} to 
all GMCs of our sample. As an example we show the results of the 
a\ plot for Orion A in Fig.[2] Similar plots for all the individual 
clouds can be found in the Appendix |Al| These plots also contain 
a cut-out of the 49 th nearest neighbour Av map, as well as the 
structure function of the cloud, determined from this map. 

The statistical error bars for the A-variance in some clouds ap- 
pear rather large. This has also been recognised by Ossenkopf et al. 
(2008a). However, they also note that (quote): "In spite of the large 
error bars, the general scaling behaviour can be accurately traced." 
Since we will measure the mass spectral index scaling exponents 
a in these diagrams we determined the weighted correlation coef- 
ficients r a (listed in the last column of Table[TJ. They show, except 



The Structure of Molecular Clouds: III 5 



Table 1. Summary of the results for the 16 investigated GMCs. The table lists the cloud name (see Paper II for coordinate ranges); the best fitting structure 
function parameters from Eq.[3]and their uncertainties (see text for details); values of the slopes s(2) and s(3); average cloud density of material with an Ay 
above 1 mag; average cloud density of material above the star formation threshold (note that 1 Mq/pc 3 corresponds to 40cm -3 ); scale of the peak in the 
A-variance analysis; scale of the peak in the A-variance analysis when applied to the star density maps (a — denotes that no peak is visible); the mass spectral 
index scaling exponent and its uncertainty from the A-variance analysis; weighted correlation coefficient for the determination of the mass scaling exponent; 
tLow density probably caused by large distance (800 pc) of cloud and hence partially unresolved higher density material. |High density due to the fact that 
the cloud is situated in a region with a general background/foreground extinction. Clouds labeled with 1 are part of the group with steep column density and 
mass distributions, and clouds labeled with 2 are part of the group with shallow distributions according to Paper II. 
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0.77 


0.07 


2.05 


0.55 


0.59 


0.13 


0.95 


0.55 


0.05 


0.01 


1.15 


1.66 
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99 




10 


0.36 


0.07 


0.980 
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0.45 


0.20 


1.18 


0.63 


0.57 


0.14 


1.82 


0.63 


0.04 


0.01 


0.70 


0.99 


12 
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0.49 


0.13 


0.944 


2 Cepheus 


0.85 


0.21 


2.31 


0.40 


0.62 


0.12 


0.69 


0.40 


0.05 


0.02 


0.86 


1.27 


5.2 


190 


16.0 




0.18 


0.10 


0.984 


1 Chamaeleon 


0.96 


0.06 


2.51 


0.34 


0.62 


0.04 


0.49 


0.34 


0.05 


0.02 


0.72 


1.05 


29 
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3.5 
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0.09 


0.949 
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0.07 
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0.38 


0.38 


0.10 
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0.38 


0.06 


0.02 


0.33 


0.44 


4.0 


42 






0.69 


0.03 


0.997 
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0.88 


0.07 


2.28 


0.45 


0.60 


0.08 


0.72 


0.45 


0.04 


0.02 


0.66 


0.93 


58 
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1.0 


1.08 


0.03 
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0.90 


0.11 
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0.47 
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0.13 
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0.47 


0.05 


0.02 


0.65 


0.92 
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100 


12.0 


4.0 


0.17 


0.03 


0.425 


2 Lupus 1 and 2 
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0.45 


0.12 


1.08 


0.47 


0.06 


0.02 


0.34 


0.47 
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680 
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0.23 


0.10 
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2 Lupus 3, 4, 5, 6 


1.08 


0.07 


2.23 


0.51 


0.49 


0.14 


0.77 


0.51 


0.06 


0.02 


-0.01 


0.02 


8.3 


1500 


3.8 




0.41 


0.04 
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0.06 
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0.45 


0.51 
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0.45 


0.05 


0.01 
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23 
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0.18 


0.04 
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1 Ophiuchus 


0.74 


0.11 


1.69 


0.54 


0.54 


0.11 


1.31 


0.54 


0.05 


0.02 


1.31 


1.89 


18 


3400 


1.5 


1.25 


0.71 


0.06 


0.975 


1 Orion A 


0.66 


0.28 


1.22 


0.55 


0.42 


0.21 


1.78 


0.55 


0.05 


0.02 


1.14 


1.58 


12 


310 


11.0 


2.0 


0.95 


0.16 


0.976 


1 Orion B 


0.83 


0.09 


2.05 


0.58 


0.57 


0.11 


0.95 


0.58 


0.05 


0.02 


1.33 


1.91 


7.0 


580 






0.52 


0.05 


0.950 


1 Perseus 


0.79 


0.12 


1.75 


0.51 


0.53 


0.08 


1.25 


0.51 


0.04 


0.01 


1.24 


1.75 


15 


1600 




1.8 


0.64 


0.12 


0.992 


1 Serpens 


0.95 


0.07 


1.63 


0.38 


0.40 


0.10 


1.37 


0.38 


0.06 


0.02 


0.17 


0.22 


100| 


860 




1.5 


0.73 


0.03 


0.957 


1 Taurus 


0.75 


0.25 


2.14 


0.55 


0.65 


0.09 


0.86 


0.55 


0.04 


0.01 


0.98 


1.41 


32 


6500 


3.5 


0.5 


0.61 


0.10 


0.950 



Table 2. Structure function parameters obtained by our test with increased 
noise. We list the noise added to the image (in units of 0.28 mag Ay, the 
1 a noise in our 49 th nearest neighbour map) and the structure function 
parameters as well as the rms. Top: Results when adding noise to our 
artificial test image. Bottom: Example for the test with a real cloud, in this 
case Corona Australis. We only used the 49 th nearest neighbour map, and 
the original image is thus represented by the first line. The noise values for 
this case hence represent the total noise, not just the added noise. 

Noise A <T/\ C ctq /3 D ar> rms oy ma 

M 

artificial test cloud 

0.0000 0.02 0.01 2.03 0.58 0.99 0.01 0.97 0.58 0.01 0.01 

1.0000 0.50 0.06 1.18 0.41 0.55 0.11 1.82 0.41 0.04 0.01 

1.0625 0.50 0.06 1.21 0.42 0.55 0.10 1.79 0.42 0.04 0.01 

1.1250 0.48 0.06 1.09 0.40 0.53 0.12 1.91 0.40 0.05 0.01 

1.1875 0.48 0.06 1.17 0.43 0.56 0.11 1.83 0.43 0.04 0.01 

1.2500 0.44 0.06 0.96 0.38 0.49 0.14 2.04 0.38 0.05 0.01 

Corona Australis 

1.0000 0.85 0.03 2.64 0.25 0.68 0.02 0.36 0.25 0.09 0.01 

1.0625 0.83 0.06 2.35 0.41 0.64 0.04 0.65 0.41 0.03 0.01 

1.1250 0.83 0.06 2.39 0.41 0.65 0.04 0.61 0.41 0.03 0.01 

1.1875 0.86 0.06 2.42 0.40 0.64 0.04 0.58 0.40 0.03 0.01 

1.2500 0.92 0.05 2.57 0.33 0.64 0.03 0.43 0.33 0.04 0.01 

for two clouds (A-Ori, Lupus 1, 2) that there is a very good correla- 
tion, with r a generally significantly larger than 0.9. 

Using our con-noise maps we investigated the positions of 
the peaks in the A-variance plots. These peaks reveal the scale at 
which a change in the structure occurs, e.g. the size of the cloud, 
the length of filaments etc.. These results are shown in Table[T] 
We find peaks in some of the clouds (Cepheus, Chamaeleon, A- 
Ori, Orion A, Lupus, Monoceros, Taurus) while for the remaining 
clouds no clear-cut peak can be identified. The identified peaks ap- 
pear over a range of scales from 1 .5 pc (Lupus 1 and 2 and Ophi- 
uchus) to 16pc (Cepheus). 



Using the star density instead of the column density maps, 
we find peaks in the A-variance for the clouds Auriga 1, Circinus, 
Corona Australis, A-Ori, Ophiuchus, Orion A, Perseus, Serpens and 
Taurus. In general, the peaks occur at smaller scales (from 0.5 pc 
in Taurus to 10 pc in Auriga 1). However, only for four clouds (A- 
Ori, Ophiuchus, Orion A, Taurus) do we detect a clear peak in both 
maps. In each case we find that the peak determined when using the 
column density map is larger, indicating that the star density maps 
trace smaller scales when the A-variance technique is applied to 
them. 

By means of the con-noise maps we also measure the mass 
spectral index scaling exponents a for each cloud. These values and 
uncertainties are listed in the last two columns of TableQ] The listed 
values are the mean of the results for the three available column 
density maps. The values for a range from 0.17 (A-Ori) to 1.08 
(Corona Australis) with a mean for all clouds of 0.54±0.27. 



3.4 The effect of noise on the A-variance 

Similar to the investigations of the structure function, we deter- 
mined the influence of noise on our analysis of the cloud structure 
with the A-variance method. The same noise increments as for the 
structure function analysis were added to the extinction maps and 
the change to the peak position of the A-variance and the mass 
spectral index scaling exponent were determined. 

Neither the peak position nor the mass spectral index seem to 
be significantly and/or systematically influenced, if an additional 
noise of the level discussed in Sect. l2.2l is added to the extinction 
maps. This applies to both, the use of our artificially generated im- 
age as well as the real cloud. Hence the obtained values can be 
considered robust. 
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4 DISCUSSION 
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Figure 2. Example of the delta variance cr^ for Orion A. We show the re- 
sults for all three available Ay maps (dash-dot 25 th nearest neighbours, 
dotted 49 th nearest neighbour, dash-dot-dot 100 th nearest neighbours). 
The corresponding vertical lines indicate the spatial resolution limit for 
each map, below which the values are meaningless. The typical scale of 
the structure in this cloud can be identified by the peak at about 1 1 pc. For 
clarity, error bars are only overplotted on the 49 th nearest neighbour data. 
The weighted correlation coefficient for the determined slope is 0.976. 



3.5 Molecular cloud density 

The determined average densities of the entire GMCs, as well as 
the densities of the material above the star formation threshold can 
be found in Table[T] 

The average densities for the entire cloud range from 
1.3 M©/pc 3 (for Monoceros) to 100 M Q /pc 3 (for Serpens). The 
low value for Monoceros might be due to small scale structures not 
being detected because of the distance to the cloud. For the Serpens 
cloud the apparently high density is due to a general offset in the 
extinction values (either foreground or background to the cloud). 
When not considering the very low and very high values, we find 
a mean value for our sample of about 15 M©/pc , which corre- 
sponds to 600 cm -3 . 

The densities of material above the star formation threshold 
are naturally higher. We find a range from 23 Mq /pc 3 (Mono- 
ceros) to 6500Mq/pc 3 (Taurus). The median value is about 
750 M Q /pc 3 or 3- 10 4 cm . Here the values for the individual 
clouds should be considered with great care, as we find larger val- 
ues for clouds with smaller distances, indicating the non-detection 
of higher column density material in more distant clouds (due to 
limited spatial resolution) and thus an underestimate of the density. 



4.1 Star Forming properties 

To investigate possible influences of cloud structure on the star 
formation properties we evaluate the current and potential future 
amount of star formation of our sample. 

i) In Paper II (e.g. bottom left panel of Fig. 4) we found that 
our sample of clouds can be divided into two groups when consider- 
ing their column density and mass distribution in the low extinction 
regions (Ay above 1 mag and below the star formation threshold). 
We indicate which cloud belongs to which group in Table[T] Clouds 
in group 1 have a shallow column density and mass distribution 
while clouds in group 2 have a steeper distribution. Hence, group 1 
members have a larger amount of material at high column densities, 
and thus more mass potentially available for star formation. One 
finds that the average fraction of mass potentially involved in star 
formation is about two to three times higher for clouds in group 1 
(see Table 4 in Paper II). 

ii) We use literature data to estimate the current number of 
young stars and mode of star formation for each cloud. In particular 
the work by Kainulainen et al. ( f2009l and Reipurth (2008a, 2008b) 
have been used. It turns out that clouds in group 1 contain typi- 
cally more than 100 young stellar objects and the majority of these 
clouds form at least one cluster of stars. Clouds in group 2 on the 
other hand, typically only form a few tens of stars and no clusters. 

Thus, group 1 contains clouds which intensely form stars, 
preferably in a clustered mode, while clouds in group 2 show 
weaker star formation activity mostly in a distributed mode. Taurus, 
without a known cluster of young stars, seems to be an exception 
as it belongs to group 1. 



4.2 Structure functions 

For our structure function analysis we selected only pixel values 
between the self shielding column density threshold of 1 mag Ay 
and the individual star formation threshold for each cloud (see Pa- 
per II), as this range represents the part of the clouds where turbu- 
lence is expected to dominate the structure. We calculated the struc- 
ture functions for each cloud using the three con-noise maps and a 
spatial scale range from 0.1 pc to 1.0 pc (traceable for all clouds). 

Using the data shown in Table[T|we investigated possible cor- 
relations between the parameters A, C against the molecular cloud 
densities (p av and psf) as well as a range of other parameters from 
Table 4 of Paper II. No correlations could be found. The analysis 
was also performed by considering two separate groups of clouds. 
There are no significant differences for any of the structure function 
parameters between the two groups. Hence, our determined struc- 
ture functions are not able to distinguish between clouds that form 
stars in clusters or in a distributed mode. 

One might expect that the values of the co-dimension C in- 
crease with distance (i.e. the fractal dimension D is decreasing), 
due to the change in appearance of the cloud. More distant clouds 
will appear to have a more simplistic structure due to details not 
being resolved. However, in our analysis no such trend could be 
found. This can be attributed to the fact that we determine the struc- 
ture functions in the scale range 0.1 pc to 1.0 pc, which is resolved 
for all investigated clouds. 

Our mean value for all clouds of A = 0.85±0. 16 is much 
higher than the values of A = expected by K41. Furthermore, it 
is inbetween the predictions of A = 2/3 (SL94, B02) and A = 1.0 
(S08). However, the scatter in our results means the value could be 
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as low as A = 0.69 or as high as A = 1 .01 . For the parameter C our 
average of 1.92±0.37 is close to the SL94 value of C = 2.0, while 
B02 theorised a value of C= 1.0. In S08 C is found to be w 1.1 for 
compressive forcing and w 1.5 for solenoidal forcing, both below 
our average value for the co-dimension. Our average is closer to the 
predicted value for solenoidal forcing (in particular considering the 
uncertainties). Hence, we conclude that the average cloud of our 
sample has a structure function hinting to solenoidal forcing of the 
turbulent field, rather than compressive driving. 

4.3 A-variance 

Using the A-variance technique, we identified peaks in plots of a\ 
against Ar to quantify the scale of structure in our investigated 
GMCs. We were only able to find peaks in a\ for about half the 
clouds. In the remainder no peaks could be identified. This could be 
caused by the fact that there are no dominant scales in those clouds, 
or that the size of the map around each cloud was too small (due 
to neighbouring clouds). Note that the A-variance will only detect 
dominant scales if they are significantly smaller than the map size. 

Scales identified average at about 5 pc. These generally corre- 
spond to the width of a filament or sheets in the cloud rather than 
the length or size of the entire cloud. In the region Lupus 3, 4, 5, 6 
the peaks correspond to the length of the different filaments. When 
we perform the A-variance analysis with the star density maps, 
dominant scales are found as well for about half the clouds. These 
scales are smaller and average at about 2 pc. This indicates that the 
star density maps trace more compact, higher extinction regions, 
generally the width of dense filaments. 

For the mass spectral index scaling exponent a, Federrath et 
al. J2010t found values of a = 0.55 for solenoidal forcing of tur- 
bulence and a = 1.34 for compressive driving. We find a consider- 
able variation in the measured values for a amongst our clouds. All 
values are below the predicted amount for pure compressive driv- 
ing. However, the mean for all clouds is a = 0.54±0.27, hence very 
close to the predictions for pure solenoidal driving of turbulence. 

When considering the average values of a for our two groups 
of clouds, significant differences can be found. While the clus- 
ter forming clouds have an average of a = 0.73±0. 18, the clouds 
with a more distributed star formation mode have a mean of 
a = 0.29±0. 1 3 (or, when excluding the clouds with low correlation 
coefficient, o? = 0.32±0.14). Based on the predictions from Fed- 
errath et al. ( 12010b . the larger value of a for the cluster forming 
clouds indicates an enhanced component of compressive driving 
compared to the clouds with a more distributed mode of star forma- 
tion. We refrain from trying to estimate a more quantitative state- 
ment about the amount of contribution from compressive driving 
in the cluster forming clouds, since the uncertainties and scatter in 
our data are clearly large. 



5 CONCLUSIONS 

We present the homogeneous investigation of a selection of 16 
nearby Giant Molecular Clouds. Utilising near infrared extinction 
maps made from 2MASS data with the near infrared colour ex- 
cess technique, we determine structure functions and perform a 
A-variance analysis. The same clouds as in Froebrich & Rowles 
(2010) were investigated. They are all reasonably nearby, away 
from the Galactic Plane, and to the best of our knowledge there 
is only one cloud along each line of sight. 

Using our extinction maps we identify two groups of clouds 



based on the slope of their column density and mass distributions of 
the turbulence dominated (low Ay ) part of the cloud. We find that 
clouds with shallower mass distributions form stars preferably in 
a clustered mode, while clouds with steep mass distributions form 
fewer stars and show a more distributed mode of star formation. 

Structure functions determined for all clouds homogeneously 
within a range from 0. 1 pc to 1 .0 pc cannot be used to distinguish 
clouds with different star formation modes. Comparing the struc- 
ture function parameters to model calculations suggests that the 
turbulence in the investigated clouds is governed preferably by 
solenoidal forcing. 

Our results of the A-variance analysis also indicate that the av- 
erage cloud in our sample is governed by solenoidal forcing. How- 
ever, clouds which form a large number of stars in clusters have an 
enhanced component of compressive driving of the turbulent field, 
in comparison to clouds with isolated star formation. Hence, com- 
pressive driving seems to lead to a more clustered mode of star 
formation. 

To quantify these qualitative findings, i.e. to determine the 
fraction of solenoidal and compressive driving in each cloud, a 
more detailed comparison of the numerical simulations with the ob- 
servational data needs to be performed. In particular viewing angle 
and resolution/distance effects in the model data need to be inves- 
tigated. 
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APPENDIX A: STRUCTURE FUNCTION RESULTS 
Al Individual clouds 
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Figure Al. Sample of the Appendix. All other figures are available online only . Top left: An extraction of the Ay map (using the nearest 49 stars) around 
Auriga 1. Gray scales are square root scaled from Omag (white) to 15 mag optical extinction (black). The size of the image at the distance of the cloud is 
133.5 pc by 70.7 pc. Bottom left: A plot showing the structure function s(p) /s(3) against p. Our data (crosses) is compared against the models of K41 (black 
solid line), SL94 (black dotted line) and B02 (black dashed line). The best fit is shown as red solid line and its parameter values (A, C, j3, D, and the rms) 
are listed in the upper left of the panel. The difference of the fit to the data against p is also shown as red solid line in the bottom of the panel. Top right: 
A-variance calculated using our con-noise maps. The dash-dot line denotes the Ay map with the nearest 25 stars used, the dotted line denotes the Ay map 
with the nearest 49 stars used, and the dash-dot-dot line denotes the Ay map with nearest 100 stars used. For clarity, error bars are only shon for the 49 stars 
data. Bottom right: A-variance calculated using our star density maps. The dash-dot, dotted, and dash-dot-dot lines are as in the upper right panel. 



